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ABSTRACT 

For  many  years  the  cylindrical  dipole  and  circular  loop  have  been 
thought  of  as  the  basic  radiating  elements.  Due  to  its  relative 
simplicity  the  study  of  the  cylindrical  dipole  is  quite  complete,  but 
owning  to  the  degree  of  complexity,  progress  in  the  theory  of  the  loop 
has  been  much  slower.  In  recent  years  the  isolated  circular  loop 
antenna  of  arbitrary  size  has  been  the  subject  of  mathematical  and 
experimental  analysis  [7],  [6].  It  was  not  until  1970  that  the  first 
complete  and  general  mathematical  solution  for  the  mutual  admittance 
of  two  arbitrary  size  coaxial  circular  loop  antennas  was  published. 
The  problem  was  formulated  by  Adler  in  terms  of  coupled  integral  equa- 
tions for  the  current  distributions  on  the  loops  [1].  The  equations 
were  solved  using  the  Fourier  series  representation  of  the  current 
distributions.  The  study  presented  here  is  the  next  logical  step 
beyond  the  work  of  Adler.  Using  his  current  equations  the  current 
distributions  are  calculated  at  each  point  around  the  two  loops.  The 
input  admittance,  electric  field  and  power  gains  are  calculated. 
Radiation  patterns  and  the  loop  current  distributions  are  plotted  for 
two  cases:  (1)  both  loops  driven  and  (2)  one  loop  acting  as  a  short 
circuited  parasite. 
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I.   INTRODUCTION 

A.   HISTORY  OF  WORK  ON  THE  LOOP  ANTENNA  PROBLEM. 

In  1897  Pocklington  excited  a  thin  loop  by  a  plane  wave  and 
obtained  a  solution  for  the  induced  current  in  the  form  of  a  Fourier 
series  expansion  [9].  Hall  en  continued  the  work  by  producing  an 
expression  for  the  current  distribution  in  the  driven  case  [3].  Due 
to  singularities  arising  in  the  Fourier  series  coefficients,  Hallen's 
technique  applied  only  to  electrically  small  loops.  Storer  attempted 
to  surmount  the  limitations  of  Hallen's  work  and  derive  an  expression 
for  the  current  in  the  case  of  a  large  loop  [10].  He  expressed  the 
series  in  integral  form  and  evaluated  the  integral  in  the  sense  of  the 
Cauchy  principal  value.  Experiments  of  Kennedy  showed  Storer' s  cal- 
culations produced  an  unacceptable  deviation  between  calculated  and 
measured  susceptance  [6].  Wu  reviewed  the  work  of  Hall  en  and  Storer 
and  discovered  that  the  shortcomings  arose  from  over-simplifying  the 
problem  [11].  That  is  they  assumed  the  electric  field  intensity  on 
the  surface  of  the  conductor  was  uniform  around  the  circumference  with  a 
filamentary  current  flowing  at  the  center  of  the  conductor. 

Wu  then  modified  Storer1 s  solution  by  changing  the  one  dimensional 
integral  equations  to  two  dimensional  integral  equations  in  order  to 
account  for  the  electric  field  variation  around  the  circumference  of 
the  conductor.  Wu  assumed  a  longitudinal  current  flow  on  the  surface 
of  the  conductor  and  allowed  no  circumferential  component.  He  showed 
this  to  be  equivalent  to  a  uniformly  distributed  current  density  on  the 
cross  section  of  the  conductor,  which  allowed  the  current  to  be  treated 


as  filamentary  and  concentrated  at  the  center  of  the  conductor.  After 
this  simplification  the  equations  were  again  one  dimensional  but  no 
longer  required  uniform  electric  field  around  the  circumference.  King 
applied  Wu's  modification  to  calculate  the  admittance  of  a  single  loop 
[7].  This  proved  Storer's  work  to  be  valid  but  only  for  the  case  of 
small  loops. 

The  first  attempt  at  solving  the  multi-loop  problem  was  by  Iizuka 
for  the  special  case  of  two  identical  loops  [5].  In  1970  Adler  advanced 
King's  work  by  devising  a  completely  general  mathematical  solution  to 
the  mutual  admittance  of  large  arbitrary  size  coaxial  loops.  Adler  also 
measured  the  admittance  in  the  two-loop  case  and  found  close  agreement 
with  theoretical  values  [1]. 

B.   SCOPE  AND  LIMITATIONS  OF  THE  STUDY 

This  paper  completes  the  study  of  large  arbitrary  size  coaxial  loop 
antenna  arrays.  Current  distributions  and  radiation  patterns  are 
calculated  and  plotted  for  various  array  geometries.  Although  only  the 
two-loop  situation  is  presented  herein,  the  theory  applies  to  the 
n-loop  situation  as  well.  Since  the  <{>  (polarized)  mode  electric  field 
component  is  of  primary  interest,  this  mode  is  examined  in  detail.  The 
6  mode  is  calculated  and  one  representative  field  pattern  is  plotted. 
The  range  of  loop  dimensions  selected  is  for  the  resonant  and  near 
resonant  case,  a,  b,  and  c  are  the  conductor  radius,  loop  radius  and 
loop  spacing  respectively;  the  range  of  parameters  are  expressed  as 

0.5  s  kb  *  2.0,     0.5  s  kc  <:  7.0 

where: 

«,   2n 


ka  is  determined  by  the  thickness  parameter: 

fl  =  2  In  2nb  and  \   is  the  wave  length, 
a 

Having  7.5  s;  ^  ^  15  represents  a  ratio  of  b/a  from  approximately  10  to 
300.  The  above  parameter  limitations  are  a  result  of  the  approximations 
used  in  the  theoretical  analysis  of  the  Fourier  coefficients  of  the 
current  distributions. 

Although  perfect  conductors  in  free  space  are  specified,  the 
effect  of  lossy  media  can  be  included  merely  by  considering  the  propaga- 
tion constant  complex.  The  loops  are  driven  by  delta  function  voltage 
generators,  symbolized  by  V  6(<f>).  The  two  coupled  antennas  are  treated 
as  a  two-port  network  and  conventional  circuit  analysis  techniques  are 
employed  to  produce  a  convenient  matrix  of  equations. 
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II,  THEORY 

The  conductor  radius  is  assumed  to  be  small  relative  to  the  loop 
radius,  which  allows  the  current  density  in  each  loop,  Ii(<fr)  and  Io(<J>)» 
to  be  approximated  by  an  equivalent  one-dimensional  filamentary  current 
distribution  at  the  center  of  the  loop.  The  electric  field  is  derived 
by  the  scalar  and  vector  potential  in  terms  of  the  distribution  of  the 
current.  Readers  are  referred  to  Adler's  work  for  the  mathematical 
details  of  the  current  equation  formulation.  Only  an  outline  of  the 
solution  is  presented  here. 

As  shown  in  Figure  1,  the  cylindrical  coordinates  for  points  on 
loops  1  and  2  are  (b-j ,  <f>',o)  and  (b2 ,4>  ,c)  respectively.  The  angle  $' 
represents  the  coordinate  of  the  current  element,  while  the  angle  $ 
represents  the  electric  field  observation  point  coordinate.  The 
distance  from  the  source  current  element  on  one  loop  to  the  electric 
field  observation  point  on  the  other  loop  is: 


R12  ■  R2]  =  V(bl~b2)2  +  c2  +  4b1b2sin2  fa^j   . 

The  distance  from  one  point  to  another  on  loop  1  is: 

Rn  =V  4b2  sin2k^\   +  4a2  sin2 

as  shown  in  Figure  2.  The  same  expression  applies  to  loop  2  if  the 
subscript  1  is  changed  to  2.  $   is  the  angular  distance  around  the  cir- 
cumference of  the  conductor  surface.  In  order  to  simply  the  notation 
in  this  development,  the  notation  for  harmonic  time  variation,  eJa)t 
has  been  suppressed. 


D  7.CK. 
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Figure  1.   Geometry  of  Two  Unequal  Loops 
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Figure  2.   Geometry  of  Points  on  a  Loop 
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e  is  the  total  electric  field  at  a  point  due  to  loop  currents 
I-.U)  and  ^U).  The  boundary  conditions  on  e  at  the  conductor 
surface  require  the  tangential  component  E^  to  be  zero  at  the  surface; 
since  no  circumferential  current  is  assumed,  the  E  component  is  the 
only  tangential  component  present.  The  delta  function  driving 
generator  placed  in  the  loop  causes  a  singularity  in  the  electric  field 
at  the  gap  such  that: 


/   E*  bl 
-n  J        *1  ' 


and  hence: 


d*  =  -v1 


Vl 


1     bl 


■* 


The  boundary  conditions  on  both  loops  are  analogous,  e.  is  the  $ 
component  of  the  total  electric  field  due  to  the  currents  on  both  loops. 
These  boundary  conditions,  when  applied  to  the  expressions  for  the  total 
electric  field  intensity,  specify  the  following  coupled  complex  integral 
equations  for  the  unknown  current  distributions: 

Y^U)  =  ^2.   J      [K11U-#')I1(*')  +  KU[W,)I2(41)W  (2) 

/n 
CK21  (*-cf> ' )  I-,  (cO » )  +  K22U-<f>')I2U')]<V       (3) 

In  these  equations 

K4i(+-+')  ■  [kbib.  cosU-*')  +  I  1-  ]  W. ,(♦-♦•)  (4) 

■J  J  K  d(^d  'J 
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where  K. .  is  the  kernel  of  the  integral  equations  and 

'  v 

[kb.b     cos(W)  +  i-i!L] 

0  K  3<f>2 

operates  on  W..,  (<J>-4>').  Wu's  contribution  to  the  solution  resulted  in 

integrating  the  term: 

e"J'kRii  ...   i    rU   e"jkRii 


get:  !!„(♦-♦')  -'-J 


to  get:  W^l*-*')  =  ~  J       r-^~ 


where  W^  (<J>-<f>')  is  the  Fourier  coefficient  of: 

e  D  n  .  When  i  f   j 

Rii 

-ikR 
W.. (♦-+')  =  £ LL  .  (6) 

Rij 

n0  is  the  characteristic  impedence  of  free  space. 

Equations  (2)  and  (3)  are  now  solved  for  I,(<j>)  and  I«(<|>)  by  the 

orthogonal  expansion  technique.  The  kernels  of  the  integral  equations 

and  the  currents  are  specified  in  terms  of  their  respective  Fourier 

series  as  follows: 

.  -jntf)' 
I.U)  =  I       In  e  (7) 

n=-a> 

K  (*-♦')  =1    aJJ  e  (8) 

n=-°° 

By  applying  equations  (7)  and  (8)  to  equations  (2)  and  (3),  inter- 
changing the  order  of  integration  and  summation  and  performing  the 
integration  over  $,  equations  (9)  and  (10)  "are  formed: 
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w 


11    I       19    9      J"^ 

jUo  I   [  ''  i'  +  J2  i2  ]  e    =  Vl5(#) 
2  p   n   n    n   n         ' 


jn<j> 


j^o  j  [-21  jl  +  a22  j2  -j     =  v  (  j 
2  n   n   n    n   n         2  xy/ 


(9) 


(10) 


After  multiplying  both  equations  by  e~J  *  and  integrating  over  <j>  from 
-n  to  n,  the  following  simultaneous  equations  are  obtained: 


$ 

12 

an 

«1 

_  1 

jnn0 

vi 

21 

22 

ii 

h 

°n 

n 

(11) 


1  2 

Equations   (11)  are  solved  for  I'   and  IR  and  substituted  back  into 

equation  (7)   to  get: 


I   U)   =  J—^cos  n<f> 
1  Jnnon 
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'1 

"n 

Vo 

22 

2 

n 

11 

12 

a 

<*„ 

n 

n 

21 

.22 

-jn<j> 


(12) 


and 


I2U) 


jnn 


£cos  n<j> 


a 


11 

n 

21 
n 


-jncj) 


(13) 


o  n 


a 


a 


11 
n 

21 
n 


.12 


.22 


The  task  now  is  to  evaluate  the  a1J's.     The  coefficients  a"""3  are  put  in 

n  n 

terms  of  the  Fourier  coefficients  <1J  of  W1-1-(<j>-<j>' )  by  letting: 

n  '  j 
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m^U-*') 


~  ij     -jn(*-*') 

)         k      e 

n=-»    n 


(14) 


where: 


/ij  =  .ij  =  1.     fU  w,.U-*')  eJnU-*')  d(j). 
n  -n      2n    y         ^ 


(15) 


-n 


The  following  difference  relationship  is  arrived  at  after  differentia- 
tion of  equation  (15): 


„1j  =  aij  =  ^jh- 
n  -n  2 


n+1         n-1 


k      Kn 


(16) 


a       =   kb.b.    <]3. 
o  1    J       I 


For  the  evaluation  of  k]^  refer  to  appendix  A, 


(17) 


The  coordinate  system  is  as  indicated  in  Figure  3.     The  expressions 
for  the  far-zone  electric  field  and  gain  are  developed  around  Adler's 
coupled  current  equations   (12)   and  (13)   rewritten  below: 


M^-s 


V     a12 
1   ao 


Vi    a12 
I     n 


Vo  a22  -  V     a22 

11      22       12     21   +&   COs(n^       ^22-y^2T 
V  ao     "ao     ao         n_1  an  an     "an     V_ 


(18) 


I2U)  =  -J 


jnn0 


*V   V2 


all   a22  _a12  a21        jfi} 
_  ao        o  O        0 


+  2\     cos(n<j>) 


all    V] 
a21    V9 


"71     22       12     21 

an'    a£      -a^   an  j 


.(19) 


Let  Ii(<J>),   the  current  in  loop  1   conductor  element  b  d<j>'   at  the  angle  <j>' 
be  as  given  in   (18)   above  and  ^U)  be  the  same  except  for  loop  2. 
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Figure  3.     Vector  potential   for  one  loop  due  to  element  of 

current  I(V)  D  d<f>'    (a)   at  arbitrary  point  and  (b) 
at  point  in  plane  of  loop;    (c)   spherical   coordinates 
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The  element  of  vector  potential  dA  at  the  point  (x,y,z)  due  to 
I,U')  is  composed  of  two  components,  dA-jx  and  dA,  where: 

dAlx  =  "dAl  sin  *'  (20) 

and 

dAly  =  +dA]  cos  (j)'.  (21) 

So  that: 

-ikR 
dAlx  =  -  H_  1^')  e-R—  sin  <f>'  bjd*'  (22) 

and 

-ikR 

dAiy =  +  4F  Ii(*,)  hr- cos  ♦'  bld*'  (23) 

where 


=  Vz2  +  s2   =  V  [rn   cos  e)2  +  s2 


R  =  Hzc   +  sc       =  If  (r  cos  q)c   +  s^  (24) 

and 

s2  =  (r  sin  e)2  +  b2  -  2brQ  sin  e  cos  ($-$*)  .  (25) 

?    2 
When    b  <<r   ,  R  -  rQ  -b  sin  e  cos  (<{>-<}> ')  . 

2    2 

The  vector  potential  from  loop  1  at  points  where  b  <<r0  is  thus  given 

by: 


_b^£jkro  /«i    ejkbi  s1necos(W)sln#.d#.  (26) 
ix     4n   ro  o*^ 

Ar  =  4^5^    f  ",(♦•)  ejkbl   sin  6  cos<*-*'>  cos  ♦'«!♦'    •   (27) 
ly         4n       r0     0»/ 

The  same  procedure  is  followed  to  get  the  vector  potential   radiated 
by  loop  2. 
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Equation  (18)  is  substituted  into  equations  (26)  and  (27).  The 
resulting  integrals  are  solved  as  follows: 


Ii  (o)  =  - 


nc 


V2  «lZ 


u   Dn 


(28) 


where  Ii    (o)   is  the  nth  term  in  the  series  for  loop  1   current  at  4>=o 

n  ^^ 

(the  feed  point)   and  %  = /u/e  where  e  is  permittivity. 


n         n       n       ^n       n 


Qi 


V     e"jkro 
4n        r« 


Ii    (1) 


B-j  =  kb-.   sin  e     . 


(29) 
(30) 
(31) 


Express  sin<{>'   and  cos<j>'    in  exponential   form  so  that 


AC    -3 

\       j2 


a  (-<a  +  Kb> 


and 


0-1 


where 


_!1  (K    +  K.) 
]v        j2         a         b 


2n  j(n+l)<j)'  jB  cos(<j>-<j>') 
Ka        /       e  e  d«j»' 


■/ 

o 


j(n-l)*'   jB^osC^-^* ) 


d<j,'    . 


(32) 


(33) 


(34) 


(35) 


These  integrals  are  of  the  form  of  the  Sommerfeld  representation  of  the 
Bessel   function: 

jBcos$  jn$ 


rn    fZn 


e        d<£>     . 


(36) 
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Letting  $  =  *-*'   gives: 

4-n     -jn*/<2Ji-*  jBcos(*-*')     jn<j>' 

■*. 


.-n     -jn*/«2Jl-*  jBcos(*-*')     j"n*' 

*«-£-•  7    e  e    "♦'  • 


(37) 


This  requires  that 


n+1  j(n+l)*    (B^ 
and 


Ka  =  2n  j    e      Jn+1     .  (38) 


n-1  j(n-l)4    (B,) 
Kb  =  2n  j    e       J^     .  (39) 

The  components  of  the  vector  potential  due  to  the  component  I-j     in 

spherical  coordinates  are: 

At"   =  (Ai   cos  *  +  AV   sin  *)  cos  e  (40) 

x„        'v« 

and 


en     xn        'yn 


A!"   =  (-Ai  sin  *  +  AV   cos  4)  .  (41) 

\      *n       'yn 

Substituting  (30),  (34),  and  (35)  into  (40)  and  (41)  gives: 

jQl 

A[   =  -=-0-  [(K  -Kb)  cos*  -j(K  +Kb)sin*]  cose  (42) 

en   t 

jQl 

=  — ^  [Kae~J*  -KbeJ*]  cos*  (43) 

2 

Af   ■  5k  [(K  +k.)  cos*  -j(Ka  -Kb)  sin*]  (44) 

*n    2 

=  -^   (Kae-j*  +Kbej*)   .  (45) 

Substitute  (38)  and  (39)  into  (43)  and  (45)  to  get: 

A[en  =  -nQlne^+n/2)[Jn+1(Bl)+Jn_/Bl)]cose        (46) 

=  -2nQ1   eJ'n(*+n/2)^    J  <B1>  cose  (47) 

n  B1       n 

and 
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<Pn       n 


>n 


V  =  "j2n  V 


(48) 
(49) 


By  summing  all  of  the  contributing  vector  potential  components  and 
using  the  recurrence  relationships: 


and 


x   n      n+1      n-1 


&>■<*>  -  j  .W  -  j  « 

n      n-1      n+1 


the  vector  potential  at  a  point  due  to  loop  1  is; 

jn(4>+IL)  /R  x 


cose 


e    -«      e, 


and 


r        °°    r  °°  Jn(<j)+n/2)       /B  v 

AVi  \  =  "j2ni  Qln  e  J"    ]       • 


Substitute  Q 


b  -jkr0 
1  =  rlu__Le H  into  equations  (50)  and  (51). 


4n2 


A'  ■  -2n  Qt  I 


AJ  =  -j2n  Q1  I 


jnU+i)    J  (kb,sine) 
2  n    n 


kb 


cose  . 


1    sine 


J'  (kb,  sine)  . 
n    1 


(50) 


(51) 


(52) 


(53) 


Since  o^  =  c^1  and  J.n^x^  =  (-l)n  ^XK   equations  (52)  and  (53)  c 
be  rewritten  as: 


an 
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AL  = 


-2n  Qt  I  n 
1 


jn(kb1  sine) 


cose 


D  kb-j  sine 


jnU+n/2)  (-l)n  -jnU+n/2) 
[e  e  ] 


(54) 


and 


Al 


=  -j2n  Q, 


* 


Vl  ao2 


J(^(kb1  sine)  +  I 
1 


12 
1   un 


V,   a 


V2«n2 
D_ 


Jr^(kb1   sine) 


jn(<fr+n/2)  _    -jn(<j>+n/2 

•  [e  +(-l)n  e  ] 


(55) 


The  electromagnetic  field  from  loop  1    is  then  given  by: 


eF    =    -ju(e  Ar    +  $  a[  )    . 


(56) 


e  ■$ 

King  shows  that  for  a  single  loop  of  moderate  size,    |kb|<2.5,   the 
significant  contributions  are  from  the  first  three  terms  of  the  series 
[7]. 

The  electric  field  patterns  for  the  <j>  and  e  modes  for  loop  1   are 
readily  obtained  using  the  first  three  terms  of  the  respective  series. 
The  identities: 

and 


,(x)  =  j   (x)  _  J,   (x) 


are  used, 
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n=0: 


k\     m  0;  AV   -  j2n  Qt 


12 


V.  a 


V2a02 
D 


J..(kb,sine) 


(57) 


£lo  =  *  El*0  =  2n  Ql  w 


Vl  *02 
V2^2 


J-j  (kb-isine) 


(58) 


(Note  that  for  the  constant  current  distribution  situation  the  series 

has  only  the  n=0  term.) 

n=l: 

Jl(kb1sin9>     cose  sin*  (59) 


A*      =  4n  Qt 


*1 


4n  Q1 


kb-jsine 


J.!    (kb,sine)  cos<(» 


(60) 


*Ef. 


1       ul 


Ef      =-j4n  Qr 


Vi   a 


12 


1   al 
Di 


Jl(kblsine)     cose  si 


kb-jsine 


n<j> 


(61) 


(62) 


E[ 


*1 


=  -j4n  (^ 


0) 


vl  «1 


12 


.22 


'2  al 


D 


n=2: 


A[      =  j4n  Ql 


v      22 
Vl   a2 

V2a22 


J'    ( kb-,   sine)   cos<j> 


2J2(kblsine)     cose  sin2<{, 
kb-.  sine 


(63) 


(64) 
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AV      =  j4n  Q1  e 


J2* 


V  a12 

VI  a2 

V2af 


J2  (kb-j  sine)  cos2<j> 


(65) 


=  e  E^  +  *  Ef, 


(66) 


:L  =  (4n  Q1U) 


vi  42 

Do 


2J2(kb1sinQ)  cose  sin2* 
kb-.  sine 


(67) 


[       =    (411  Q,a)) 
4>2 


v,42 

u       22 

V2  a2 

Do 


J£  ( kb-,   sine)  cos2<|) 


(68) 


The  field  from  loop  2  is  found  in  the  same  manner.  To  find  the 
total  field  at  a  point  for  each  mode,  each  term  of  the  series  for  loop  2 
for  both  modes  is  multiplied  by  the  phase  factor,  exp  (jkc  cose)  to 
account  for  the  separation  of  the  loops.  The  total  electric  fields  for 
the  <{>  mode  and  e  mode  are  expressed  as  follows: 


and 


e£  =  E\    +   Ej   e 
6    %  \ 


jkc  cose 


+  E  .  +  E£  e 


jkc  cose 


<J>-|   s^ 


jkc  cose 


+  EtT   +  Eo   e 

M>2    <f>2 


jkc  cose   r    r   jkc  cose 
4  =  Ele  E2e  e      +Ele,+  E2eie 

°o   eo  1    °1 


_r 

'2q2 


>e~  £( 


jkc  cose 
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The  gain  of  the  array  is  calculated  by  dividing  the  square  of  the  field 
by  the  power  radiated.  Where: 


Power  Radiated  =  Prad  =  V*  ■ 


IN 


V?  •  Y 


IN' 


(69) 


In  order  to  find  the  input  admittance,  it  is  necessary  to  first  derive 
expressions  for  self  and  mutual  admittance.  By  treating  the  two  coupled 
loops  like  the  two-port  network  in  Figure  4,  equations  (70)  and  (71)  can 
be  written  directly: 

1,(0)  =  Yn  V,  +  Y12  V2  (70) 


I2(o)  =  Y21  V,  +  Y22  V2  . 


(71) 


s 

A 

Figure  4.  Two-Port  Network 


Ii(o)  and  Io(o)  are  the  feedpoint  currents.  Because  the  network  is 
reciprocal,  Y^  =  Y2-, .  Solving  for  mutual  admittance: 


Mo) 


12 


V1=0      J^o 


Self-admittance  is  by  definition: 


11 


22 


V0) 

Io(0) 


v2=o 


V^O   . 


v2=i 


(72) 


(73) 


(74) 


24 


The  input  admittance  is  by  definition: 

h  v2 

Y  =  -J1    =   Y        +  _£.     Y 

\  Yll  Y12 


(75) 


Y1N     "         =  Y22  +  —    Y21      * 
2       V2         "       V1       *' 


(76) 


In  order  to  simplify  the  digital   computer  program  used  later  in  solving 
for  the  electric  field  and  gain,   it  will   be  assumed  that  both  driving 
voltages  are  purely  real.     This  will   allow  Yt«  to  be  replaced  by  G^ 
where  G™  equals  the  real   part  of  YI[Nj. 
Then: 


'IN 


v2 

"  Gll   +  —  G12 
Vl 


and 


GINo  =  G22  +  „     G21 


(77) 


(78) 


Finally, 


Gain,   = 


and: 


Gain, 


*     +E? 


"1* 


rad 


r  r 

El     +E2 


rad 


(79) 


(80) 


r '  r '  r  r 

E£     and  E2     represent  E2  and  E2     respectively  multiplied  by  the  phase 

+jkc  cose 

factor  (e  )  which  accounts  for  effects  of  loop  separation. 
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III.  DISCUSSION  OF  NUMERICAL  TECHNIQUES  AND  COMPUTER  PROGRAM 

The  a*3   terms  were  computed  on  an  IBM  System/360  Model  67  digital 
computer  in  FORTRAN  IV  language  using  Adler's  basic  program  [1].  Beyond 
this  point  the  programming  was  straight  forward.  Flow  charts  on 
Figures  5a,  5b,  and  5c  show  the  major  steps  in  the  program.  The  flow 
chart  of  Figures  5a  and  5b  is  for  computing  I-jU),  IoU)>  Y-,,  Ys2>  Y-j-j, 
Y22>  and  Y^  where  Y,  and  Y  «  are  the  se^  admittance  of  loops  1  and  2 
respectively  with  mutual  coupling  effects  ignored.  Figures  5a  and  5c 
show  the  steps  in  computing  Gain^  and  Gaine.  The  computer  results  for 
current  in  the  isolated  loop  were  compared  with  King's  [8],  while  the 
gain  was  compared  with  the  work  of  Harrington  and  Mautz  [4].  In  both 
cases  the  agreement  was  exact.  The  incrementing  of  parameters  b  and  c 
in  the  program  allows  for  changing  loop  size  and  loop  spacing. 
Incrementing  e  and  <j>  allow  for  computation  of  gain  patterns. 
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INPUT: 
bj  ,  b^ ,  c  , 

Ab  ,  A.  c  , 


CALCULATE 

II    12    22 

<       »  K   »  K~ 
o     o     o 


K   =  NEW 
n 

k   ,=  OLD 
n-  I 

SAVE=  NEW 


n  =  n+1 


< 


SAVE  = 


<n 


n-1 


=  K 


a 


CALCULATE 

12    22 


a 


a 


CALCULATE 
11    12 


22 


NEW  =  K 
n 

n  =.  n  ~  1 


a 


K   =  SAVE 
n 

CALCULATE 

11    12    22 


a. 


a 


n  '  "n  '  "  n 
OLD   =  SAVE 


Figure  5a.  Flow  Chart  of  Computer  Program 
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CALCULATE 

JsVjs2»YlT 

Y22'Y12 
n=n+l 


$ 


'OUTPUT  LOOP 
DIMENSIONS 

^si'Ys2'3irY22>/ 

Y12 


CALCULATE: 


$ 


XL 


OUTPUT 

Isl(*),Is2(*)»I1(*).I2(f) 

m=l ,2, . .  ,m 


Figure  5b.  Flow  Chart  of  Computer  Program 
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4 


b2  =  initial  h^ 
c  =  c  +  AC 


STOP 


$ 


Figure  5b.  Flow  Chart  of  Computer  Program  (Continued) 
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0 


SET  e 


CALCULATE 
GAIN^ 

GAINfl 


OUT 


P! 


UT 


GAIN.}) 
GAINQ 


<f)    =    <j>    +    Acf) 


<D 


<D 


^0 


e  =  e  +  Af 


Figure  5c.  Flow  Chart  of  Computer  Program 
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4 


bp  =  b?  +  Abp 


b2=INITIAL  b2 


C  =  C  +  AC 


* 


STOP 


Figure  5c.  Flow  Chart  of  Computer  Program  (Continued) 
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IV.  DISCUSSION  OF  RESULTS 

A.   CURRENT  DISTRIBUTION  AND  GAIN  PATTERN  CORRELATION 

Figures  6  and  7  show  the  patterns  of  Gain,  in  the  x-z  plane  and  the 
x-y  plane  respectively.  In  all  cases  presented  kb-j=l  and  loop  1  is  in 
the  x-y  plane.  Loop  2  is  positioned  coaxial ly  above  loop  1  at  a 
distance  c.  Figures  8,  9,  10,  11  and  12  show  the  corresponding  current 
distributions  for  the  case  of  both  loops  driven  in  phase  by  one  volt 
delta  function  generators. 

For  .5  <  kc  <  1  and  .5  <  kb~  <  1.5  the  radiation  patterns  do  not 
change  significantly  and  only  a  slight  change  in  phase  of  loop  2  current 
is  noted.  When  kc  >  1  radical  changes  in  pattern  shapes  occur.  This 
shows  that  as  the  separation  becomes  larger  than  about  X/6,  the  effect 
of  the  array  factor  predominates.  When  kb2  =  .5  the  circumference  of 
loop  2  is  A/2.  In  this  case  the  loop  supports  a  1/2  cycle  standing 
current  wave. 

Figures  13  and  14  are  the  Gain,  patterns  in  the  x-z  and  x-y  planes 
respectively,  for  the  case  when  loop  2  is  parasitic.  Figures  15,  15, 
17,  and  18  are  the  corresponding  plots  of  the  current  distributions. 
When  the  two  loops  are  closely  coupled,  kc  =  .5  and  resonant  the 
currents  are  large  and  almost  equal  in  magnitude  but  180  degrees  out  of 
phase.  Even  though  Gain<|>  is  relatively  high  in  this  configuration,  the 
benefit  of  large  currents  is  not  realized  because  of  the  phase  differ- 
ence between  the  two  currents.  As  the  separation  of  the  loops  is 
increased  from  approximately  A/13  to  approximately  A/6,  the  currents 
drop  by  a  factor  of  about  five;  however,  the  Gain  drops  only  slightly. 
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kc 


.5 


6.3 


kb. 


.5 


1.5 


Figure  6.     Gain  $  patterns   in  x-z  plane  ($=  o°,180°);   feed  point 

at  <j)=o°,  e=90°;  <j>=o°  at  right,   e=o°  at  top;  n=10,   kb1=l, 
Vi=V2=l;   kc=.5,l,3,6.3;   kb2  =    .5,1,1.5,2. 
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kc 

kb2 
.5 


.5 


6.3 


1.5 


Figure  7.  Gain  <j>  patterns  in  X-Y  plane,  e=90°;  feed  point  at 
bottom  (<j>=o°);  ft=10,  kb-,=l,  Vi=V2=l;  kc  =  .5,1,3, 
6.3;  kb2  =  .5,1,1.5,2. 
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kc 
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6.3 


kb, 


.5 
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Figure  13.     Gain  $  patterns   in  X-Z  plane  (^=0° ,180°) ;  feed  point 
at  (j)=o°,   6=90°,   4>=o°  at  right,   8=o°  at  top;  n=10, 
kb-|=l ,  V^l,  V2=0;   kc=.5,l,3,6.3;   kb2=.5,l  ,1 .5,2. 
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Figure  14.     Gain  <j>  patterns   in  X-Y  plane,   e=90°;   feed 
point  at  bottom  (<j>=o°;  Q=10,   V -j  =  1 ,  V2=o, 
kb-j =1 ;   kc=  .5,1 ,3,6.3;   kbo^.S.l ,1 .5,2. 
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This  observation  is  explained  as  follows:  although  the  currents  are  much 
smaller  when  c  -   x/6,  the  phases  of  the  two  currents  are  no  longer 
exactly  opposing.  This  in  turn  allows  the  total  vector  potential  of  the 
<J>  mode  in  both  cases  to  be  almost  equal. 

B.   CURRENT  DISTRIBUTIONS 

Referring  again  to  the  current  plots,  the  following  observations 
are  made.  When  both  loops  are  resonant  and  driven,  the  current  peak 
increases  as  c  increases  from  X/12  to  x  and  decreases  as  c  increases 
beyond  x.  When  c  is  increased  beyond  2x  (not  shown)  the  mutual  coupling 
is  negligible  and  the  loop  currents  are  very   close  to  those  of  isolated 
loops.  This  phenomenon  is  due  to  the  behavior  of  the  input  admittance, 
where: 


YIN1  "  Yll  +  Y12 


Vl  =  V2=  1- 


By  computing  Y-,-,,  Y^o*  Y-io  and  yin-i  9   tne  following  was  noted:  (a)  with 
close  coupling  Y,2  has  a  large  negative  real  part  which  approaches  zero 
as  c  is  increased;  (b)  the  real  part  of  Y-j-j  is  always  positive;  (c) 
Y j ^  increases  as  c  increases.  Therefore,  since  1(0)  =  Y^V,  the  current 
increases  as  Yj^  increases. 

Of  particular  interest  is  the  behavior  of  the  current  in  the 
parasitic  loop.  In  all  cases  both  resonant  and  non-resonant  the 
magnitude  of  the  parasitic  current,  although  reduced,  paralleled  the 
magnitude  of  the  current  on  the  driven  resonant  loop.  In  the  resonant 
case  when  the  spacing  was  close,  the  two  currents  were  180  degrees  out 
of  phase,  but  as  the  spacing  approached  x,  the  currents  approached  the 
in-phase  condition.  Also  in  the  parasitic  case  the  phases  of  both 
currents  tend  to  be  constant  with  an  abrupt  reversal  at  the  nulls. 
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C.   GAIN  PATTERNS 

Figures  6  and  7,  the  Gairu  patterns  in  the  x-z  plane  and  x-y  plane 
respectively  are  again  referred  to.  In  these  figures  ft=10,  V-j=V2=l  volt, 
kb-|=l,  kc=.5,  1,  3  and  6.3  from  left  to  right,  while  kb2=.5,  1,  1.5  and 
2  from  top  to  bottom.  In  all  patterns  presented  "+"  marks  the  antenna 
center  and  M- — -"  marks  unity  gain.  With  kb2  fixed,  the  patterns  are 
seen  to  change  shape  as  kc  increases  from  .5  to  6.3.  This  is  due  in 
small  measure  to  the  change  in  mutual  coupling  and  primarily  to  the 
array  factor.  The  smallest  Gain,  in  the  z  direction  occurs  at  kb-j  = 
kb2  =  1  and  kc  =  3  (c  -  X/2);  the  largest  Gain,  in  the  z  direction 
occurs  also  with  the  resonant  array  but  with  kc  =  6.3  (c  -  x) .     In  the 
former  case  the  vector  potentials  were  nearly  180  degrees  out  of  phase 
and  in  the  latter  both  vector  potentials  were  in  phase.  In  the  extreme 
case  where  kb2  =  2  and  kc  =  6.3,  lobing  occurs  due  to  the  combination  of 
the  array  factor  and  the  diverse  current  distributions  on  the  two  loops. 

Gain^  for  the  parasitic  case  is  illustrated  by  Figures  13  and  14. 
The  higher  gains  and  more  interesting  patterns  are  attained  in  this 
configuration.  As  would  be  expected,  less  lobing  now  occurs  and  the 
parasitic  element  functions  as  a  reflector.  Of  special  interest  are  the 
kb2  =1.5,  kc  =  1  and  kb2  =  2,  kc  =  6.3  patterns.  In  the  first  of 
these  particular  arrangements  the  loop  array  pattern  resembles  that  of  a 
log  periodic  antenna.  Figure  21  shows  the  effect  of  changing  kb2  to  1.2 
and  varying  kc  from  .6  to  1.2.  When  kb2  =  kc  =  1.2  the  antenna  is  uni- 
directional. In  contrast  to  this  configuration,  note  the  resulting 
pattern  when  kb2'  =  2  and  kc  =  6.3.  The  pattern  is  almost  square  with 
the  gain  across  the  top  nearly  uniform.  Figure  14  shows  a  section  of 
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these  same  patterns  in  the  X-Y  plane.  These  patterns  closely  resemble 
those  of  the  familiar  single-loop  pattern  in  shape. 

The  effect  of  changing  the  b/a  ratio  is  illustrated  in  Figure  19. 
The  change  in  G  is  from  7.5  to  15,  which  represents  a  b/a  ratio  change 
of  about  10  to  300.  Although  the  effect  is  measurable,  it  is  concluded 
that  the  loop  array  is  not  sensitive  to  changes  in  ft. 

Figure  20  shows  the  consequence  of  ignoring  the  mutual  coupling 
when  calculating  the  gain  pattern  (the  dotted  curves  are  with  coupling 
ignored).  In  Figure  20a  (kc  =  .5  and  kb-,  =  kb^  =  1)  the  coupling  is 
very  significant,  but  as  seen  in  Figures  20b  and  20c  (kc  =  6.3,  kb2  =  1 
and  1.5  respectively)  the  coupling  has  little  effect. 

Sample  patterns  of  Gaine  are  shown  in  Figure  22.  Although 
significant,  this  mode  is  certainly  of  secondary  importance. 
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kb  .5  1  1.5  2 

2 


Q 


7.5 


15 


Figure  19.     Gain  <j>  patterns  in  X-Z  plane  (<j>=o°,180°)   for 
V1=V2=kb1=l,  kc=6.3;   kb?=.5,l ,1 .5,2;  fi=7.5,15; 
<{i=o0  at  right,   e=o°  at  top. 
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kc 
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1.5 
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1 


Figure  20.  Gain  <j>  patterns  in  X-Z  plane  with  (solid  lines) 
and  without  (dotted  lines)  mutual  coupling 
effects.  n=10,  kb-j  =  1 ,  <j>=o°  at  right,  e=o°  at 
top. 
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kc 
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1.2 


Figure  21.     Gain  4>  patterns  in  X-Z  plane  (<J>=o0,180°)   for 

=10,  V^l,  V2=o,   kb-,=l ,   kb2=1.2,   kc=.6,.8,l  .2, 


¥2 
kb, 

kc 
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1 
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1.5 
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(a) 
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Figure  22.  Gain  e  patterns  for  fi=10,  kb,=l  in  (a)  <f>=45°, 
225°  plane;  (b)  <j>=90°,  270°;  (c)  <j>=90\  270°. 
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V.  SUGGESTIONS  FOR  FURTHER  STUDY 

Since  this  study  is  of  a  theoretical  nature  an  immediate  require- 
ment is  to  confirm  the  results  by  experiment.  First  the  current  distri- 
butions and  next  the  field  intensity  patterns  for  various  arrangements 
should  be  measured.  Also  of  interest  in  the  two  loop  array  is  the 
consequence  of  open  circuiting  the  parasitic  loop. 

The  current  and  gain  equations  then  should  be  solved  for  the  three 
element  case.  Of  interest  would  be  the  consequence  of  arranging  three 
or  more  elements  to  make  a  log  periodic  loop  antenna. 
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APPENDIX  A 

Evaluation  of  Fourier  Coefficients  <1J 

n 

Case  I:   i  =  j 

Combining  equations   (5)  and  (15)   results   in: 


ii           !         AII/-II     J'ne     -J'kRii 
K      =  '     '       e        e .      .  (81) 


WJ 


^)ZJJ  ~R7—  d*  d« 


n 


By  substituting  A.  =  2aj  sin^-  in  equation  (81)  and  interchanging  the 
order  of  integration  one  arrives  at  the  expression: 

ii    7      /*2ai     W 
«"  =  £   /   ]  1_L_  dA,  (82) 

n  J  r— 1 

M     -A2 


where: 


1    r11  jne  -k  R^ 


Mn(A.)  =  —  /  e    e "  de;  n  >  o  (83) 

nv  i;   2n07        r.. 

n 

Storer  showed  that  for  the  case  of  A.<<b.: 


l   /8bi\  ?   n"1 
lni    n   \Ai  /   n  ^ 


2m+l 
i  /   "  m=i 

i    r2kb] 

2      I      [n2n(x)  +  j  J2n(x)  ]  dx  (84) 

o  J 

where  B  (x)  is  the  mth  order  Lommel -Weber  function  and  Jm(x)  is  the  mth 
order  cylindrical  Bessel  function  of  the  first  kind. 
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Wu  devised  an  approximation  method  to  evaluate  equation  (84)  based 
on  the  facts  that  (a)  the  integral  representation  of: 
2    „  1 


I  K0{Ai(^-^J 


where  K  (x)  is  the  modified  Bessel  function  of  the  second  kind  of  zero 
order  and  the  integral  representation  of  b^M  (A-)  are  similar  and  (b)  for 
A.j  small  the  difference  between  b-Mn(A-)  and  the  small  argument  approxi- 
mation for  KQ  is  independent  of  A^ .  The  consequence  of  this  is  for 

moderate  size  n, 

.2,2 


1        h  _  *h 


<VF 


wv-^  2KolAiV?" k  I  +lnV"^ 


2  Cn  -n    /   [C2n(x)  +  j  J2n(x)  ]  dx  (85) 


n-1  , 
where  C  »  y  -  2    £      +  ln(4n)  and  y  is  Euler's  constant. 


m=o  2m+1 


Because  A-  is  small,  the  K  term  in  equation  (85)  can  be  slightly 
modified  with  a  negligible  loss  in  accuracy  to: 


By  employing  equations  (85)  and  (86)  equation  (82)  can  be  integrated 


to  get: 


ii  „   1   (  2  ,  r^A   I  (^1\   +  2  C „  -         (87) 


Kn     "   2nb.  )  u   ,xo  \b- 

r2kbi 

n  J     [fi2n(x)  +  j  J2n^x^  dx 
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and: 


Ko1  =  2lbT    \]n(ir)    -no/^o(x)  =  J  J2n(x)3  dxj    .  (88) 

Case  II:   i7j 

By  combining  equation  (15)  with: 

Rid  =  [^ibj  sin2  (iiil)  +  (bi   -  bd)2  +  c2]2 

and 

-jkR-- 

Kij 
equation  (89)   follows: 

n  9        o   o  1 

fc1J  _  1_  f    eJ"VJk  [^b,.  sin^Mbj-bj)2^2^  rf 

11   2n_/    [4bibj  sin2i+  (brbj)2  +  c2]l       *' 

It  is  necessary  to  solve  equation  (89)  by  numerical  methods  and  the 
reader  is  referred  to  Adler's  original  work  for  the  details  of  this 
process  [1], 
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13   ABSTRACT 

For  many  years  the  cylindrical  dipole  and  circular  loop  have  been  thought  of 
as  the  basic  radiating  elements.  Due  to  its  relative  simplicity  the  study  of  the 
cylindrical  dipole  is  quite  complete,  but  owning  to  the  degree  of  complexity,  pro- 
gress in  the  theory  of  the  loop  has  been  much  slower.  In  recent  years  the  isolated 
circular  loop  antenna  of  arbitrary  size  has  been  the  subject  of  mathematical  and 
experimental  analysis  [7],  [6].  It  was  not  until  1970  that  the  first  complete  and 
general  mathematical  solution  for  the  mutual  admittance  of  two  arbitrary  size 
coaxial  circular  loop  antennas  was  published.  The  problem  was  formulated  by  Adler 
in  terms  of  coupled  integral  equations  for  the  current  distributions  on  tne  loops 
[1].  The  equations  were  solved  using  the  Fourier  series  representation  of  the 
current  distributions.  Tne  study  presented  here  is  the  next  logical  step  beyond 
the  work  of  Adler.  Using  his  current  equations  the  current  distributions  are 
calculated.  Radiation  patterns  and  the  loop  current  distributions  are  plotted  for 
two  cases:  (1)  both  loops  driven  and  (2)  one  loop  acting  as  a  short  circuited 
parasite. 
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